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1  Overview 


Project  title:  Radar  Control  —  Optimal  Resource  Allocation 

This  report  addresses  several  issues  arising  in  discussions  with  Braharn  Hirned  as  well  as  our 
work  to  date. 

The  key  components  are: 

1.  Track  before  detect  techniques:  This  arose  in  a  discussion  with  Braharn  Himed.  The 
description  given  here  is  for  information  only.  It  is  a  technique  we  have  implemented 
several  times  in  the  past,  alongside  other  innovations,  as  have  several  other  researchers. 
This  part  of  the  document  has  been  presented  in  previous  reports  to  AFRL  as  part  of  a 
Prometheus  contract. 

2.  Mode  scheduling:  again  this  is  a  presentation  of  earlier  work  by  the  group,  also  reported 
under  a  Prometheus  contract,  and  again  in  response  to  a  request  by  Braharn  Himed. 

3.  Opportunistic  bandwidth  radar:  this  is  new  work  aimed  at  attempting  to  understand 
scheduling  in  a  time-dependent  bandwidth  limited  situation.  The  idea  is  to  choose  pieces 
of  bandwidth  as  and  when  they  become  available. 

4.  Clutter  modelling  for  detection:  this  represents  part  of  our  work  on  the  contract.  A 
component  of  waveform  and  other  resource  scheduling  requires  clutter  to  be  modelled,  in 
real  time.  Here  is  presented  a  simple  approach  that  appears  to  perform  better  than  some 
existing  techniques. 

5.  Information  gain  from  radar  scheduling:  this  is  an  early  attempt  to  understand  how  to 
compare  gain  in  a  radar  from  measuring  immediately  rather  than  waiting  to  choose  a 
better  waveform.  While  this  work  is  preliminary  it  suggests  that  it  is  not  always  better  to 
schedule  resources  for  gain  of  information. 

6.  More  simulation  based  studies  were  carried  out  and  results  are  reported  in  Section  7. 


2  Integrated  detection  and  tracking 

Tracking  algorithms  in  radar  typically  work  with  a  set  of  “detections”  obtained  after  applying 
statistical  tests  to  the  radar  returns  in  each  cell  of  the  surveillance  area.  The  statistical  testing 
procedure  involves  comparing  the  matched  filter  output  to  a  threshold  selected  to  give  a  specified 
false  alarm  probability.  This  traditional  approach  does  not  consider  the  dependence  of  tracking 
performance  on  threshold  selection.  This  dependence  was  first  considered  in  [1]  where  the 
detection  threshold  was  selected  to  minimise  an  approximation  to  the  mean  square  error.  A 
similar  idea  was  pursued  in  [2,  3].  All  of  these  methods  are  based  on  the  sub-optimal  probabilistic 
data  association  filter  (PDAF)  [4]  and  the  analysis  performed  in  [1].  The  basic  idea  is  as  follows. 
The  PDAF  provides  a  Gaussian  approximation  to  the  posterior.  Let  P^\t  denote  the  covariance 
matrix  of  this  approximation  for  the  state  at  time  k  based  on  measurements  up  to  time  t.  It 
was  shown  in  [1]  that,  averaged  over  the  measurements,  the  one-step  ahead  covariance  matrix 
can  be  approximated  as 

E[-Pt|t]  ~  Pt\t-i  ~q(T,S)Wt  (1) 

where  q(T,S)  G  [0, 1]  is  the  information  reduction  factor  for  a  detection  threshold  T  and  SNR  S 
and  Wt  is  a  positive  definite  matrix.  Assuming  the  SNR  S  is  known,  the  trace  of  the  expected 
covariance  matrix  can  be  minimised  by  maximising  q(T,  S)  as  a  function  of  T.  Although  q(-,  S) 
cannot  be  evaluated  in  closed-form  it  was  shown  in  [2]  that,  for  a  special  case,  a  function  can  be 
fitted  to  q( •,  S )  and  the  optimal  value  of  T  found  in  closed- form.  Interestingly,  the  optimal  T  for 
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tracking  purposes  usually  results  in  a  false  alarm  probability  far  in  excess  of  those  commonly 
used  in  radar  signal  processing. 

A  different  approach  was  used  in  [5]  where  the  available  knowledge  of  the  target’s  position 
was  used  to  formulate  a  Bayes  detector  which  increases  the  detection  probability  in  the  vicinity 
of  the  estimated  target  position.  This  results  in  a  detector  with  a  spatially-varying  threshold. 

The  existing  methods  described  in  [1,  2,  3,  5]  have  a  number  of  limitations: 

•  All  methods  are  based  on  the  sub-optimal  PDAF.  It  would  be  of  interest  to  consider  the 
problem  of  integrated  detection  and  tracking  using  better  alternatives  such  as  particle 
filters. 

•  The  methods  of  [1,  2,  3]  are  based  on  MSE  minimisation.  This  is  not  always  of  interest 
when  tracking  in  clutter.  In  such  cases,  track  continuity  may  be  of  more  interest  than 
tracking  accuracy. 

•  All  methods  are  limited  to  the  case  where  it  is  known  that  only  a  single  target  is  present.  In 
general,  the  number  of  targets  cannot  be  assumed  to  be  known  and  methods  for  handling 
multiple  targets  are  required.  Extending  existing  methods  to  this  case  is  of  great  interest. 


3  Mode  Scheduling  and  Simulation 

This  section  for  the  most  part  is  an  abbreviated  version  of  old  reports,  submitted  as  part  of  a 
Prometheus  contract,  we  make  no  intellectual  property  claims,  but  the  ideas  here  appear  to  be 
relevant  to  a  discussion  between  Bill  Moran  and  Braham  Hirned. 

3.1  Mode-Scheduling  for  Radar 

Modern  radar  system  has  some  specific  objectives: 

1.  detection  of  all  targets  in  the  region  being  illuminated,  i.e  surveillance 

2.  estimation  of  the  position  and  trajectory  of  those  targets,  i.e  tracking 

3.  identification  of  the  targets, 

4.  assessment  of  threat  posed  by  each  of  the  targets,  etc. 

Accomplishing  each  of  those  aims  depends  on  the  particular  mode  that  the  radar  system  is 
using.  We  envisage  a  radar  in  which  there  is  some  significant  control  over  the  radar  system  at 
a  fast  enough  rate  that  it  is  possible  to  change  its  mode  of  operation  on  a  time  scale  consistent 
with  objective  tasks  above.  Moreover  we  assume  that  these  changes  will  effect  the  performance 
of  the  radar  in  accomplishing  this  task.  We  assume,  in  order  to  be  specific,  that  the  task  is 
estimation,  which  might,  as  in  the  context  of  the  simulator  we  have  built,  be  a  “track-before- 
detect”  capability  and  so  encompass  both  Obj-1  and  Obj-2.  However,  we  will  describe  the 
objective  in  the  language  of  tracking.  We  will  also  assume  that  the  mode  being  controlled  is 
one  or  more  of  any  number  of  radar  parameters,  such  as  Waveform,  PRI,  CPI,  Beam  shape, 
Beam  direction,  Frequency,  etc.  The  problem  now  becomes  a  control  problem:  how  to  optimally 
schedule  the  control  parameter  so  as  to  maximize  a  cost  function  which  measures  the  effectiveness 
of  the  radar  in  accomplishing  its  task. 

Typically,  an  MoE  will  have  some  estimate  of  the  current  clutter  situation  as  well,  at  least 
in  a  tracking  or  track-before-detect  situation,  of  the  targets.  Thus  it  is  possible  to  simulate  the 
effects  of  any  given  control  (such  as  waveform)  on  the  scene  and  judge  its  capability  according 
to  a  cost  function.  This  sets  the  scene  for  consideration  of  the  problem  as  one  in  control  - 
in  fact  stochastic  control.  We  note  that  classical  control  techniques  assume  a  knowledge  of 
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the  state  of  the  system.  Here  we  only  have  an  estimate  of  the  state  of  the  system  based  on 
the  prior  measurements.  This  estimate  is  often  referred  to  as  the  information  state  and  can 
often  be  considered  as  a  replacement  for  the  state  in  making  a  stochastic  control  problem  into 
a  deterministic  one.  There  are  results  (so-called  separation  theorems )  that  show  that  under 
a  wide  range  of  circumstances  the  information  that  arises  from  prior  measurements  can  be 
condensed  down  to  a  “sufficient  statistic”  that  carries  all  of  the  information  needed  from  the 
prior  measurements  needed  by  the  scheduler  in  arriving  at  its  optimal  schedule. 

3.1.1  One-step  ahead  control 

There  are  many  references  in  the  literature  to  adaptive  radar  techniques.  Almost  all  of  these 
describe  either  passive  adaptive  radar  methods  (in  which  the  adaptation  is  done  only  on  receive 
in  the  processing  of  the  return)  or  active  adaptive  methods  in  which  the  scheduling  is  done  in 
a  one-step  ahead  or  greedy  way.  This  method  seeks  at  the  next  PRI  (or  whatever  the  time-scale 
used)  to  choose  the  waveform  that  is  optimal.  This  may  seem  to  be  the  obvious  thing  to  do. 
One-step  ahead  scheduling  will  certainly  buy  more  than  fixed  mode  radars  but  multi-step  ahead 
scheduling  will  buy  more  and  there  are  circumstances  in  which  what  is  paid  in  computation  is 
worth  it. 

3.1.2  Closed  Loop  Scheduling 

This  is  the  optimal  multi-step  scheduling  technique.  At  each  step  the  cost  function  is  evaluated 
for  the  next  N  steps  based  on  the  predictions  described  above.  The  optimal  schedule  is  selected 
and  implemented  for  the  next  epoch  only ,  and  the  process  is  repeated  using  the  new  measurement 
to  enhance  the  prediction.  In  this  way  all  information  is  used  at  each  step  of  the  scheduling 
algorithm. 

As  can  be  imagined,  this  is  computationally  very  demanding.  The  calculation  of  the  optimal 
schedule  is  a  difficult  dynamic  programming  problem  and  this  has  to  be  done  at  each  epoch  (so 
perhaps  each  PRI)  and  in  order  to  do  this  the  new  measurement  has  to  be  incorporated  into 
the  estimates  and  predictions  at  each  epoch. 

3.2  Computer  Simulations 

Under  the  Prometheus  effort  we  developed  a  computer  GUI  that  allowed  simulation  of  various 
tracking  scenarios,  i.e  the  number  of  targets  on  the  scene  with  various  birth/death  times,  their 
trajectories,  radar  cross-sections,  etc.  Various  clutter  scenarios  are  provided,  ex.  stationary  of 
moving  clutter,  land  or  sea  clutter  and  the  amount  of  clutter  on  the  scene  could  be  selected. 
Waveform  (using  a  small  library  of  waveforms)  and  beam  scheduling  algorithms  are  implemented 
to  allow  longer  time  for  performing  the  surveillance  (smart  surveillance).  The  simulated  results 
are  shown  in  the  GUI  and  the  statistic  of  the  outcome  is  computed  and  shown.  An  example  of 
such  outcome  is  shown  in  Figure  1.  Here,  the  three  scheduling  algorithms  are  compared,  one 
and  two  steps  ahead  scheduling  with  smart  surveillance  were  tested  against  the  algorithm  with 
surveillance  scheduled  at  regular  time  intervals. 

4  Opportunistic  Bandwidth  Radar 

4.1  Background 

In  the  situation  when  the  whole  bandwidth  is  unavailable,  the  available  smaller  bandwidths  can 
be  utilized  with  opportunistic  radar.  Suppose,  that  the  total  bandwidth  is  1  GHz,  but  only  IMhz 
of  it  is  available  at  any  time.  In  this  work  we  assume  that  available  bandwidths  are  detectable 
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1  Step  Ahead 

2  Steps  Ahead 

No  Waveform  Scheduling 


Average  RMSE  26.177  m;  21.4915  m;  23.5269  m; 


Figure  1:  The  RMSE  and  Target  Revisit  Time  with  Smart  vs.  Regular  Surveillance 


and  known,  but  appear  at  random  frequencies  within  the  total  frequency  band.  Either  up-sweep 
of  down-sweep  LFM  waveform  can  be  transmitted  in  the  available  band. 

s  =  [exp(2vri^^/2),  exp(2vr  A^/2),  . . . ,  exp(27ri^/2)],  (2) 

where  n  is  the  number  of  pulses  transmitted,  B  is  either  1  GHz  or  —  1GHz  total  band  for  up-sweep 
or  down-sweep  LFMs,  L  -  total  pulse  length.  t\,  t2,  ■  ■  ■  ,tn  are  computed  as  hiliJi+Ml ;  where 
/i,  ■  •  • ,  fn  are  frequency  bands  available  at  time  instances  j  =  1 , .n  and  A /  =  1MHz  available 
bandwidth.  Square  brackets  represent  vector  concatenation,  curly  brackets  represent  (frequency) 
interval.  A  spectrogram  for  such  transmission  is  shown  in  Figure  2.  The  autocorrelation  of  such 
signal  is  shown  in  Figure  3.  It  can  be  seen  that  it  has  good  autocorrelation  properties.  Hence,  a 
stationary  scatter  could  be  detected  with  such  a  signal.  This  is  not  the  case  for  moving  scatterer, 
because  the  return  pulse  will  have  an  unknown  doppler  shift.  This  is  shown  in  Figure  4. 

4.2  Signal  processing 

It  is  possible  to  improve  the  performance  by  rearranging  the  returns.  For  stationary  reflectors, 
the  return  signal  can  be  “unscrambled”,  by  sorting  fi,  f2,  ■  •  ■ ,  fn  in  acceding  order  and  concate¬ 
nating  the  returned  signals  in  this  order.  The  auto-correlation  of  this  signal  is  shown  in  Figure  5. 
If  the  Doppler  shift  is  known,  then  the  above  procedure  can  be  performed  ,  after  “Doppler 
compensation”,  i.e  by  modulating  the  return  with  exp(— 2-Kifdt). 

We  consider  the  following  algorithm. 

•  Perform  standard  Doppler  processing 

•  Compensate  return  signal  with  Doppler,  obtained  in  previous  step 

•  Unscramble  the  signal 
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Figure  2:  A  spectrogram  of  the  waveform  transmitted  with  an  opportunistic  radar 
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Figure  4:  A  processed  return  of  the  waveform  transmitted  with  an  opportunistic  radar  reflected 
from  a  moving  scatterer 
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Figure  5:  An  auto-correlation  of  the  waveform  transmitted  with  an  opportunistic  radar,  un¬ 
scrambled 
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•  Detect  Range 

•  Use  this  information  to  track  the  targets 

•  Refine  the  search  from  step  1,  using  the  optimum  waveform,  computed  on  the  tracking 
step. 

The  simulation  results  of  the  above  algorithm  is  shown  in  Figure  6.  In  this  example  target  speed 
10.5  m/s.  This  shows  that  the  proposed  method  significantly  reduces  the  sidelobes  of  the  return 
even  in  the  presence  of  the  Doppler  shift. 
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Figure  6:  Processed  Radar  Return  with  and  without  Doppler  Compensation  Method  with  slow 
moving  target  (target  speed  is  10.5  m/sec) 


5  Sea  clutter  modeling 

Sea  clutter  is  the  unwanted  echoes  resultant  from  radar  returns  reflected  off  the  ocean  surface 
when  operating  a  radar  within  a  maritime  environment.  The  prevalence  of  this  physical  phe¬ 
nomena  routinely  presents  a  challenge  for  the  radar  designer,  frequently  presenting  itself  as  the 
dominant  noise  component  in  particular  at  low-grazing  angles.  It  has  the  potential  to  place 
limitations  on  the  detectability  of  targets  that  may  be  obscured  its  presence  sea  clutter  due  to 
the  physical  structure  of  the  sea. 

We  are  interested  in  developing  parsimonious  models  which  exploit  the  inherent  structure 
of  sea  clutter.  When  the  sea  is  considered  to  be  fully  developed,  sea  clutter  over  a  reasonable 
observation  period  (>  30  s)  can  be  observed  to  possess  a  quasi-periodic  structure  for  a  given 
patch  of  ocean.  This  phenomenon  is  synonymous  with  the  observed  waves  present  at  the  ocean 
surface,  with  a  characteristic  periodic  structure  clearly  visible  in  a  radar  range-time  plot.  By 
assuming  that  the  sea  clutter  is  a  cyclostationary  process,  we  estimate  a  set  deterministic  com¬ 
ponents  that  form  the  texture  or  macrostructure  component  of  the  sea  clutter  which  has  been 
shown  to  valid  model  for  longer  observation  periods.  Previous  statistical  estimates  of  these  de¬ 
terministic  components  [6,  7]  did  not  consider  the  spatial  relationship  of  adjacent  rangebins  in 
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the  measurement  set,  instead  independently  estimating  the  texture  for  each  range  bin.  Here  we 
seek  to  derive  both  the  temporal  and  spatial  relationship  of  the  estimated  texture  components. 

5.1  The  Compound  Gaussian  Sea  Clutter  Model 

At  any  point  of  the  ocean  surface,  the  observed  waves  are  the  result  of  the  locally  generated  wind 
(capillary)  waves  (short  waves)  modulated  onto  the  longer  (gravity)  waves  that  may  last  several 
days  and  propagate  long  distances.  This  permits  the  representation  of  the  sea  clutter  using  two 
different  surface  roughness  scales:  one  that  describe  the  small  scale  roughness  or  speckle  often 
associated  with  whitecaps  and  foam  and  the  large  scale  roughness,  the  texture  associated  with 
the  swell.  Then,  the  complex  envelope  of  the  return  can  be  modelled  as,  for  n  =  1, . . . ,  N, 

z(nT)  =  \J~Tn  xn  (3) 

where  T  represents  the  radar  pulse-repetition  interval  (PRI)  that  is  fixed  over  the  interrogation 
period,  {rn}  is  the  texture  and  {xn}  is  the  speckle.  Here  we  further  assume  that  thermal  noise 
is  negligible  in  comparison  to  the  clutter.  In  the  compound  Gaussian  model,  the  speckle  {xn} 
is  a  circular  complex  white  Gaussian  process. 

The  conventional,  non-parametric  way  of  estimating  the  sea  clutter  is  to  apply  a  low  pass 
filter.  Let 

7(nT)  =  \z{nT)  |2  =  rn\xn\2  (4) 

Then,  we  can  estimate  the  texture  non-parametrically  using  a  moving  window  (MW)  filter  [8] 
of  length  L  +  1  as 


n+L/2 

Tn  =  1/(L  +  1)  ^2  1  (mT) 

m=n—L/2 
n+L/2 

=  Tnj  (L  +  1)  ^2  \xn\2  =  t na\  (5) 

m=n—L/ 2 

assuming  that  thermal  noise  is  negligible.  In  order  for  the  approximation  of  (5)  to  be  valid,  we 
select  L  such  that  L  \  <C  L  <C  LT,  where  Lx/T  and  Lr/T  are  the  speckle  and  texture  correlation 
periods,  respectively  For  X-band  sea  clutter,  coherence  analysis  shows  a  decorrelation  period  of 
about  5  ms  to  15  ms  for  the  speckle,  whereas  the  texture  decoration  times  are  in  the  order  of 
several  seconds  [9,  10]. 

5.2  Parametric  texture  estimation 

A  parametric  estimator  of  texture  in  the  compound  Gaussian  model  (3)  can  provide  a  more  suit¬ 
able  representation  of  sea  clutter  for  the  purposes  of  adaptive  sensing.  The  observed  periodicity 
of  waves  suggests  a  parameteric  texture  model  of  the  form  [11]: 

K 

Tn  (Ok)  =  +  +  22  Ak  cos  (27rfknT  +  <f>k)  (6) 

k= l 

where  Ok  =  {A,/,  (/>}*  for  the  I\-  th  order  spectral  components  of  the  texture.  We  further 
assume  that  {+}j[L0  are  deterministic  positive  constants,  {/+(+  G  (0,cj/27t]  and  {+}(L0  G 
(— 7T,  vr]. 

Previously  the  model  (6)  has  been  considered  independently  for  each  range  cell.  However,  in 
general,  we  expect  that  for  a  fully  developed  sea,  the  wavelengths  of  the  cosine  components  of 
the  swell  are  similar  and  the  waves  travels  ‘together’  in  a  characteristic  periodic  behaviour.  This 
spatial  relationship  can  be  exploited  by  assuming  the  existence  of  a  linear  phase  relationship  in 
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range  between  the  phases  of  each  frequency  component.  Consider  a  texture  model  that  is  now 
the  planar  sum  of  a  set  of  harmonics  across  of  a  series  of  contiguous,  evenly  spaced  and  sized  R 
range  bins. 

K 

Tn,r  (Ok)  =  A0  +  ^2  Ak  cos  ( 2nfknT  +  4>k  +  ipkr)  (7) 

k= 1 

where  now  6  k  =  {A,  /,  (/>,  ip}o  are  the  common  harmonic  components  of  the  texture  across  all 
of  the  range  cells  and  {riPk}k=o  ^  (— 7r,  7t]  is  used  to  describe  the  phase  shift  of  the  k- th  harmonic 
of  the  texture  in  range.  This  variable  is  expected  to  be  also  be  related  to  the  wavelength  of  the 
harmonic  component  itself. 

This  finally  allows  us  to  now  describe  the  sea  clutter  as  a  function  of  both  range  and  time: 

Zr(nT)  —  sj Trir  J'ri.r  (8) 

where  {xn;r}  are  independent,  zero-mean  circular  complex  Gaussian  random  variables. 

Texture  estimation  from  the  measured  intensity  requires  the  retrieval  of  closely  spaced  har¬ 
monic  components  in  the  presence  of  the  speckle,  a  complex  multiplicative  noise.  The  closeness 
of  these  components  prevents  us  from  using  traditional  spectral  analysis  methods  to  estimate  the 
harmonics.  Instead  we  use  a  modified  version  of  the  CM-RELAX  algorithm  [6]  that  recursively 
estimates  harmonics  of  the  cyclic  mean  from  the  texture  intensity.  This  RELAXation  algorithm 
was  originally  introduced  as  an  optimisation  approach  for  estimating  harmonic  parameters  in 
data  sequences  corrupted  by  autoregressive  (AR)  noise  [12,  13,  14]. 

We  consider  two  variants  of  the  CM-RELAX  Newton.  The  first  independently  estimates 
the  frequency  components  of  each  range  bin,  as  in  [6].  Since  we  wish  to  establish  a  spatial 
relationship,  the  frequency  estimates  from  each  range  bin  are  then  clustered  to  determine  the 
spatial  frequency.  The  second  method  estimates  the  frequency  components  for  R  adjacent  range 
bins  simultaneously  by  the  exploiting  the  spatial-temporal  relationship  of  the  sea  clutter,  as 
in  (7).  We  refer  to  the  first  algorithm  as  independent  CM-RELAX  and  the  second  as  joint 
CM-RELAX. 

The  method  we  introduce  here  differs  from  CM-RELAX  in  a  number  of  ways.  Firstly,  we 
exploit  the  symmetry  of  the  intensity  periodogram  when  estimating  its  frequency  components 
through  maximisation.  The  prevents  the  need  for  estimate  both  sides  of  the  periodogram  to 
be  separately,  significantly  reducing  the  time  complexity  required  for  the  algorithm  to  run. 
Additionally,  we  use  the  Newton-Raphson  method  to  find  the  maximum  of  the  periodogram  of 
the  k- th  order  of  the  cyclic  mean.  This  numerical  method  provides  a  practical  means  of  finding 
closely  spaced  frequency  components  without  needing  to  perform  extremely  long  FFTs.  By 
doing  this,  we  found  that  it  was  unnecessary  to  continuously  re-estimate  the  already  estimated 
frequency  components  until  the  texture  estimate  settled  as  previously  done  in  CM-RELAX. 

5.3  Validation 

The  sea  clutter  data  used  to  validate  the  algorithm  was  collected  at  Osbourne  Heads  Gunnery 
Range,  Dartmouth,  Nova  Scotia,  Canada  by  the  McMaster  University  Intelligent  PIXel  (IPIX) 
X-band  Polarimetric  Coherent  Radar  during  the  OHGR  -  Dartmouth  Campaign  as  shown  in 
Figure  7.  This  radar  is  an  experimental  frequency  agile,  high-resolution  originally  developed 
for  the  purpose  of  iceberg  detection  research.  It  features  a  8kW  TWT  operating  at  8.9  GHz  to 

9.4  GHz.  The  parabolic  dish  antenna  has  a  beamwidth  of  0.9°  with  a  cross-polarisation  isolation 
of  greater  than  33  dB. 

The  radar  was  located  on  a  cliff  overlooking  the  Atlantic  Ocean  at  a  height  of  30  m  in  a 
staring  configuration  above  mean  sea  level  at  an  azimuthal  bearing  of  129.5°  over  the  open 
ocean. 

We  selected  a  measurement  performed  on  the  November  11th  1993  at  16:36:25,  corresponding 
to  the  data  file  19931111_163625_starea.  cdf .  At  the  time  of  collection,  the  sea  was  reported 
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Figure  7:  Radar  data  collection  site  at  OHGR,  Nova  Scotia.  Source:  [15] 


to  having  a  significant  wave  height  of  0.7  m  with  a  average  wave  period  of  4.4  s.  During  the 
measurement  period,  the  radar  was  operating  at  9.39  GHz  with  a  PRF  of  1000  Hz  with  a  range 
resolution  30  m,  sampled  at  15  m  over  a  period  of  131  s.  We  utilise  range  bins  1  to  6.  This 
equates  to  ranges  of  2574  m  to  2649  m  at  a  resolution  of  30  m,  sampled  at  15  m  due  to  the 
placement  of  a  target  at  2660  m.  The  results  reported  here  correspond  to  measurements  made 
in  a  VV  polarisation. 

The  results  shown  here  apply  to  the  CM-RELAX  Newton  algorithm  being  run  over  a  a 
data  set  of  length  120  s.  For  the  intensity  data  Yk-\  FFT  computation,  we  found  that  an  FFT 
length  of  218  and  28  provided  sufficient  resolution  for  a  good  initial  peak  estimate  of  u>k  and  c fa , 
respectively.  For  comparative  purposes  with  the  non-parametric  MW-Filter  texture  estimator, 
we  have  used  a  length  of  210  samples  with  equates  to  a  period  of  1.024  s. 

We  first  consider  the  correspondence  between  the  spatio-temporal  frequencies  obtained  by 
the  independent  and  joint  CM-RELAX  algorithms.  Figure  8  shows  the  frequencies  obtained 
by  independent  CM-RELAX  for  K  =  7  across  the  series  of  six  range  bins.  Five  clusters  of 
closely-spaced  frequencies  are  clearly  evident.  This  suggests  that  the  spatio-temporal  relation¬ 
ship  proposed  in  (7)  has  merit.  Two  of  the  remaining  frequencies  were  not  common  across  the 
range  bins.  This  can  be  attributed  to  the  algorithm  not  estimating  the  same  frequencies  for  all 
range  bins. 

The  phases  of  the  five  clustered  frequencies  are  plotted  as  a  function  of  range  in  Figure  9(a). 
The  gravity  wave  phase  of  the  frequency  components  found  by  joint  CM-RELAX  are  shown  in 
Figure  9(b).  A  close  correspondence  between  the  slopes  of  the  phases  in  Figures  9(a)  and  9(b) 
can  be  seen.  This  is  confirmed  in  Figure  10  which  plots  the  phase  rates,  or  spatial  frequencies, 
against  the  temporal  frequencies  for  independent  and  joint  CM-RELAX.  The  reasonably  close 
agreement  between  the  two  algorithms  agrees  with  the  notion  that  the  texture  that  we  observe 
in  any  range  bin  is  the  result  of  the  sea  swell  that  is  travelling  across  the  surface  of  the  ocean  in 
range  over  time,  as  predicted  by  (7). 

We  now  consider  the  prediction  performance  of  the  independent  and  joint  CM-RELAX 
texture  estimation  algorithms.  Performance  is  measured  by  the  RMS  error  between  the  predicted 
texture  estimate  and  the  MW  estimate,  averaged  over  time  and  range  bin.  The  results  are 
shown  in  Figure  11  for  times  up  to  10  s.  The  algorithms  perform  similarly  in  terms  of  prediction 
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Figure  8:  Estimated  harmonics  {Ok}  for  I\  =  7 
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Figure  9:  Clustered  harmonic  estimates  {Ok}r  for  K  =  7 
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Figure  10:  Derived  Phase  rate  compared  with  jointly  estimated  phase  rate  'ipx 

performance.  The  RMS  error  increases  rapidly  up  to  2  s  and  then  remains  approximately 
constant.  However,  their  behaviour  as  a  function  of  K  differs  substantially.  Joint  CM-RELAX 
performs  best  for  K  =  3  while  independent  CM-RELAX  favours  values  of  K  between  4  and  8  at 
different  times.  We  attribute  this  to  the  ability  of  the  joint  CM-RELAX  to  more  easily  identify 
the  harmonic  components  of  significance  from  the  intensity  data. 

5.4  Conclusion 

We  have  developed  a  compound  Gaussian  model  for  sea  clutter  in  which  the  texture  component 
is  modelled  by  a  sum  of  sinusoids.  In  contrast  to  previous  work  with  this  model,  we  considered 
the  spatio-temporal  frequency  of  the  texture  components.  This  is  motivated  by  the  idea  that 
waves  travel  across  space  in  a  reasonably  orderly  manner.  Comparisons  with  the  usual  method 
of  independently  estimating  frequencies  in  each  range  bin  using  real  data  demonstrated  the 
validity  of  our  model. 

The  proposed  joint,  spatio-temporal  sea  clutter  model  has  several  advantages.  First,  it  allows 
the  sea  clutter  to  be  represented  by  a  much  smaller  number  of  parameters,  therefore  resulting  in 
a  more  parsimonious  model  which  is  faster  to  work  with.  Second,  because  inference  is  performed 
in  2D  frequency  space,  there  is  greater  separation  between  the  frequency  components.  As  a  result 
the  CM-RELAX  algorithm  performs  more  reliably.  This  is  manifested  in  a  greatly  reduced  run 
time  and  the  achievement  of  good  performance  with  a  smaller  number  of  frequency  components. 
A  paper  [16]  reporting  this  work  has  been  submitted  and  is  currently  under  review. 

6  Information  Gain  from  Radar  Scheduling 

6.1  Background 

Considerable  research  effort  has  gone  into  the  algorithms  for  scheduling  of  radar  resources;  typ¬ 
ically,  waveforms,  but  also  beam  direction  and  shape  and  other  tunable  parameters  of  radars 
[17,  18].  Such  radar  resource  scheduling  usually  demands  massive  computation.  Even  myopic 
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(b)  Joint  CM-RELAX 

Figure  11:  RMS  prediction  error  for  (a)  independent  and  (b)  joint  CM-RELAX 
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scheduling  requires  evaluation  of  each  feasible  value  of  the  tunable  parameters  against  a  current 
representation  of  the  state  of  knowledge.  Even  though  parsimonious  representations  of  the  scene 
are  possible,  the  computation  is  non-trivial.  When  continuous  versions  of  the  parameters  are 
used,  some  kind  of  gradient  descent  (or  worse  still  stochastic  search)  is  deployed  also  requiring 
many  computer  clock  cycles.  Parsimony  of  representation,  while  effective  in  reducing  computa¬ 
tion  also  leads  to  reduced  validity  of  the  optimal  choice  of  radar  resource.  In  the  non- myopic 
context,  the  computational  problem  becomes  exponentially  more  difficult.  Deployment  of  an 
optimal  scheduling  policy  over  more  than  just  a  few  radar  pulses  is  computationally  infeasible 
even  on  normal  research  timescales,  and  certainly  not  at  the  rate  required  for  real-time  deploy¬ 
ment  in  pulse-Doppler  radar  system,  where  such  computation  would  require  to  be  completed  in 
at  most  a  few  microseconds. 

Leaving  aside,  for  a  moment,  the  issue  of  computer  time  needed  to  find  the  optimal  resource, 
the  same  kind  of  issue  arises  when  observing  and  tracking  multiple  radar  targets.  Time  spent 
on  a  particular  target  reduces  the  time  available  to  observe  others  and  to  detect  new  targets. 
The  same  kind  of  issues  are  at  work  in  this  context. 

Returning  to  the  computation  problem,  on  the  one  hand  then,  a  radar  system  might  wait 
to  complete  the  calculation  of  an  optimal  waveform  for  a  particular  scene,  during  which  time 
the  target  could  have  significantly  changed  its  state,  and  on  the  other  choose  to  use  its  current 
best  “guess”  with  the  limitations  of  the  use  of  a  “less  than  optimal”  radar  resource.  This  report 
explores  how  a  radar  tracking  system  can  benefit  from  spending  more  time  in  scheduling  of 
radar  resources  —  and  thereby  achieving  better  performance  —  versus  using  a  less  than  optimal 
choice  of  parameters  and  obtaining  more  timely  information. 

We  analyze  a  very  simple  situation  where  the  measurement  quality  improves  with  time  spent 
scheduling  the  measurement.  We  consider  the  situations  when  measurement  error  is  reduced 
by  scheduling  (eg.  scheduling  waveforms,  number  of  PRI,  etc),  and  when  probability  of  target 
detection  is  increased  by  scheduling  (eg.  beam  scheduling,  RM  coding,  etc).  We  model  target 
dynamics  as  Gaussian  processes  and  use  the  Kalman  filter  as  tracking  algorithm.  Covariance 
of  target  state  is  used  as  a  measure  of  improvement.  Notice  that,  in  the  Gaussian  case,  this 
measure  is  directly  related  to  information  gain. 

When  the  state  variable  is  one  dimensional  the  optimal  time  allocated  for  scheduling  can 
be  computed  analytically.  While  this  case  is  extremely  simplistic,  it  guides  our  understanding 
of  the  more  general  case.  An  analytic  solution  is  not  available  when  the  state  is  represented 
by  higher  dimensional  state  variables,  but  we  are  able  to  give  a  theoretical  bound  to  the  time 
allocated  for  scheduling. 

In  this  work  we  assume  that  the  measurement  and  track  processing  times  are  negligible 
compared  to  the  scheduling  time.  This  assumption  is  not  unreasonable  for  most  modern  radar 
systems.  Hence  we  use  the  term  “scheduling  time”  to  mean  “time  between  measurements”. 


6.2  One  dimensional  case,  Gaussian  model 

This  very  simple  case,  while  unrealistic  in  practice,  gives  some  insights  that  are  useful  for  the 
higher  dimensional  situation.  We  assume  that  the  state  variable  x  and  the  measurement  y  are 
scalars  (for  example:  in  problem  of  estimating  the  range  of  target) 

x  =  dynamic  noise 

(9) 

y  =  x  +  measurement  noise. 

The  process  is  tracked  using  a  Kalman  Filter,  so  that  after  the  nth  measurement  the  state 
posterior  estimate  is  given  by  a  Gaussian  distribution  with  mean  xn\n  and  variance  &n\n. 

When  no  measurement  is  available,  the  state  evolves  as 


*£n+l|n  xn\ni 
&n+l\n  &n\n  T  Qt, 


(10) 
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where  t  is  time  and  q  is  a  dynamic  noise  variance  for  unit  time.  That  is,  our  uncertainty  about 
the  process  is  assumed  to  increase  linearly  over  time. 

We  define  radar  regime  as  operating  in  a  “no  scheduling”  mode  when  the  measurements  are 
taken  at  the  regular  time  interval  t  =  1  epoch.  Now  let  t  be  the  spent  scheduling  before  the 
next  measurement  yn+i  is  taken. 


6.2.1  Measurement  accuracy  is  improved  by  scheduling 

Radar  range  variance  relates  to  measurement  SNR  [19]  as  r  =  ro  +  where  vq  is  a  fixed 

error  variance  (due  to  radar  receiver  noise  and  propagation  conditions)  and  Ax  is  a  size  of  range 
resolution  cell.  Consider  an  example  when  SNR  depends  of  scheduling  time  as  ct,  for  some 
positive  c,  then 

Ar2 

r  =  r„  +  ^.  (11) 

Notice,  that  c  is  SNR  at  time  t  =  1  (i.e  in  ‘no  scheduling”  mode). 

For  this  and  similar  examples,  we  want  to  write  this  expression  in  more  abstract  terms  as 
r  =  a  +  b/t,  where  a,  b  are  real  positive  constants. 

The  posterior  state  variance  is  computed  by  Kalman  filter  as 


_  (tfn| n  +  qt)(b/t  +  a) 

<-rn+l|ra+l  j  / ,  |  .  ,  ,  vT^J 

b/t  +  a  +  crn|n  +  qt 

We  are  interested  in  finding  the  scheduling  time  t  >  1  which  minimizes  <rn_|_i|n+i.  Lets  assume 
that  the  tracking  system  is  stable,  i.e  is  does  not  diverge  in  the  “no  scheduling”  mode.  The 
steady  state  variance  a  in  this  case  is 


<J  =  \fq(^a  +  4 b  +  q)/ 2  -  q/2. 


(13) 


a  is  a  lowest  possible  state  variance  in  “no  scheduling”  mode,  therefore  it  is  appropriate  to  allow 
<rn|n  =  a.  The  question  is,  at  what  time  t  >  1  posterior  state  variance  surpasses  a.  Simplification 
shows  that  <rn+1|n+1  is  a  rational  function  of  f 

_  ab  +  (aa  +  bq)t  +  aqt2  , 

an+1|n+1  “  b+(a  +  a)t  +  qt 2  '  (  } 

Denoting  numerator  by  A{t)  and  denominator  by  B(t),  we  observe  that  in  the  domain  of  interest 
(i.e  for  t  >  1)  both  polynomials  are  positive,  hence  o‘n_)_1|n+1  has  no  zeros  or  poles.  The  derivative 
of  <7n+1|n+1  is  proportional  to 

— n+^n+1  ot  A'(t)B(t )  —  B'(t)A(t)  =  q(a2  —  bq)t 2  +  2 bq(a  —  a)t.  +  b(—a2  +  bq).  (15) 

It  now  can  be  shown,  that  for  t  >  1 

•  If  a  >  \fbq ,  then  dcTn+^ri+1  is  positive,  hence  crn+1|r)+1  is  increasing,  therefore  we  lose 
accuracy  by  scheduling. 

•  If  a  =  \/bq,  then  dCTn~[j1.|ri+1  =  0,  hence  we  gain  nothing  by  scheduling. 

•  If  a  <  \/bq,  then  d<Tn+^n+1  <  hence  <7n+1|n_|_i(i)  is  a  strictly  decreasing  with  t  with 
minimum  value  a  at  infinity.  Therefore  we  improve  accuracy  by  schediding  and  scheduling 
time  is  unlimited.  Incidentally,  in  this  case  lim  crn+ i|n+i  =  a. 

Obviously,  similar  argument  can  be  applied  to  a  wide  range  of  r,  for  instance  any  decreasing 
r(t)  <  a  +  b/t  and  a  <  \fbq  will  will  also  allow  unlimited  “scheduling”  time.  Of  course,  in 
practice,  the  scheduling  time  will  have  to  be  restricted  by  the  target  range,  target  detection, 
presence  of  clutter,  etc.  These  would  have  to  be  taken  into  consideration  when  designing  a  real 
radar  tracking  system. 
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6.2.2  Probability  of  detection  is  improved  by  scheduling 

Now,  we  will  analyze  the  scenario,  when  the  scheduling  improves  probability  of  target  detection 
Pd ■  The  posterior  variance  is  computed  as 


_  °n+l|rc 

®n+ 1  n+1  ^n+l|re\-^  Pd  .  “ 

<-rn+l|n  T  1 


(16) 


Consider  an  example  when  r  is  a  constant  and  pd  =  1  —  a/t,  where  0  <  a  <  1. 

As  before  we  use  steady  state  variance  for  “no  scheduling”  as  starting  point.  In  this  example 


is 


yfq  (q  +  4?’  —  Aar)  —  q  +  2  a  q 
2  —  2  a 


(17) 


We  want  to  find  a  scheduling  time  t*  >  1,  which  minimizes  Equation  (16)  and  provides 
<7n+1|n+1  <  a.  The  derivative  of  crn+1in+1  is  proportional  to 


da . 


n+l  |  n+l 

dt 


oc  q(r2  +  aq(r  —  a))t 2  —  2 aa2qt  —  aa2(a  +  r). 


(18) 


This  is  quadratic  polynomial  with  at  least  one  real  positive  root,  hence  it  has  the  solution  t* 
which  is  computed  as 


a  a2  q  +  a  r  y/aq  (a  +  r  +  aq) 
aq2  r  —  1  a  a  q2  +  q  r2 


(19) 


We  are  only  interested  in  t*  >  1  and  solving  the  inequality  we  obtain  the  following. 


•  If  a  <  0.5  and  r  <  (2 a-i )2  ’  then  t*  >  1,  and  the  highest  accuracy  is  achieved  at  at  t*. 

•  If  a  >  0.5  and  r  >  then  t*  >  1,  and  the  highest  accuracy  is  achieved  at  t*. 

•  Otherwise  -  longer  scheduling  time  always  results  in  loss  of  accuracy. 


Examples  of  the  state  variance  as  a  function  of  scheduling  time  for  the  two  cases  are  shown 
in  Figures  12  and  13. 


6.2.3  Probability  of  detection  is  improved  by  scheduling,  another  example 

Consider  another  example  of  the  above  problem.  Lets  suppose  that  at  time  t  =  1,  pd  =  a  <  1. 
Lets  also  suppose  that  after  time  T  of  scheduling  the  probability  of  detection  is  pd  =  1.  We 
would  like  to  find  an  answer  to  a  question:  what  is  maximum  value  of  a  for  which  there  will  be 
an  improvement  in  state  variance  at  time  T.  We  have  t  =  T 


(a  +  Tq)r 
<7T'T  -  a  +  r  +  T, 


(20) 


where  a  is  the  steady  state  variance  in  ‘no  scheduling”  mode  (i.e  at  t  =  1) 


a  = 


q  —  2  aq  +  (q(q  +  Aar))1/2 
2  a 


The  solution  is  the  following 


(21) 


(2  qr  —  3Tq2  +  2  Tr2  —  qu—  2  ru  +  2  q2  +  T2q2  —  AT qr  +  2  T2qr  +  Tquy/Tq  +  2TruyfTq) 

2 (Tq  -  q  +  Tr)2 


(22) 

where  u  =  y/Tq  (4r  +  Tq).  For  example,  for  q  =  1,  r  =  10,  T  =  10  the  probability  of  detection 
in  “no  scheduling”  mode  a  has  to  be  smaller  than  0.3332.  For  q  =  1,  r  =  10,  T  =  3  a  <  0.5659. 
Results  for  other  values  of  T  are  shown  in  Figure  14. 
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Variance  of  the  Posterior  State  Distribution,  case  :  measurement  noise  variance  r=a+b/t 


Figure  12:  State  variance  as  a  function  of  scheduling  time  when  measurement  accuracy  is  im¬ 
proved  by  scheduling.  In  this  example  q  =  1,  a  =  1,  b  =  10 


Variance  of  the  Posterior  State  Distribution,  case  :  Probability  of  detection  pd=1  — c*/ 1 


Figure  13:  State  variance  as  a  function  of  scheduling  time  when  probability  of  detection  is 
improved  by  scheduling.  In  this  example  q  =  1,  a  =  0.7,  r  =  10 
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Figure  14:  Maximum  pd  at  t  =  1  to  permit  improvement  by  scheduling  at  time  t  =  T.  In  this 
example  q  =  1,  r  =  10 


6.3  Two  Dimensional  Case 


Suppose  now,  that  the  range  and  Doppler  are  measured,  the  state  dynamics  and  measurements 
are  described  as 


x  =  (/  —  F)  x  +  dynamic  noise 
y  =  x  +  measurement  noise, 


(23) 


where  /  is  the  identity  matrix,  F 


\ 

I  and  r  is  time.  The  dynamic  noise  is  a  zero  mean 


^  ,  the  measurement  noise  is  drawn  from 

a  Gaussian  pdf  with  zero  mean  and  covariance  R.  In  the  first  example  we  assume  that  the 
pd  =  1.  Let  R(t )  be  a  measurement  noise  covariance  matrix  which  can  improve  with  scheduling 
time  t.  Unlike  in  one  dimensional  case,  the  improvement  could  be  achieved  in  size,  shape  and 
orientation  of  R.  The  aim  is  to  minimize  the  determinant  of  the  posterior  state  covariance 
matrix  as  a  function  of  scheduling  time.  We  write 


Gaussian  process  with  covariance  Q(t')  =  q 


Pn+l\n  ~  PPn\n,PT  +  Q(t) 

det(Pn_|_i|n+1)  det(I  T>ri_|_^|n(i:>n_|_^|ri  -F  R(t ))  )  det(-Pn_|_i|n)) 


(24) 


where  t  is  the  scheduling  time.  As  above  it  is  (almost)  equal  to  the  time  between  nth  and 
(n  +  l)th  measurements. 


6.3.1  Measurement  Accuracy  improves  with  Scheduling  Time 

In  this  example  we  postulate  that  scheduling  will  improve  measurement  accuracy  as  R  =  RO/t , 
then  the  calculations  show,  that 


det(Pn+1|n+1)  —  |Po|  > 


(25) 
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where  A^(t)  =  aot4  +  ait3  +  a2t2  +  a^t  +  a±  and  B$(t)  =  b^t6  +  bit 5  +  b2t4  +  b^t3  +  b^t2  +  b^t  +  b$ 
are  positive  (for  t  >  1)  4th  and  6th  order  polynomials  in  t,  with  following  coefficients 

[q2AP22q,  I2pi2<?,  I2png,  12|P|]  (26) 


for  A4  and 


q2, 4p22q,  4q(3p12  +  r22),  12(pnq  +  p22r22  ~  qr12), 

12(|P|  +  2r22Pi2  -  2p22r12  +  qrn),  12tr(PPy 1)|-R0|,  12|P0| 


(27) 


for  Bq. 


det(Pn+1|n+1)  is  positive  on  the  domain  of  interest  (f  >  1)  and  it  has  no  zeros  or  poles  in  the 
domain.  It  behaves  like  1/f2  for  large  t.  So  we  conclude  that  there  are  no  local  minima  in  the 
domain. 

Numerically  simulations  also  show  that  det(Pn+i|n_|_i)  is  decreasing  with  time.  In  this  exam¬ 


ple  q  =  1,  R  = 


5  2 
2  1 


It  appears  that  for  reasonable  values  of  parameters  det(Pn+1|n+1)  is  a 


Figure  15:  Determinant  of  State  Covariance  as  a  function  of  scheduling  time  when  measurement 
accuracy  is  improved  by  scheduling.  In  this  example 

strictly  decreasing  function  of  time,  hence  is  no  limit  on  maximum  allowed  time  for  scheduling. 

6.3.2  Probability  of  detection  is  improved  by  scheduling 

In  this  scenario  we  assume  R  to  be  a  constant  with  time  and 

Pd  =  1  —  a/t,  0  <  a  <  1.  (28) 

We  have 

|Pn+1,n+1|  =det((I-pdK)(FPFT  +  Q))  =  (29) 

In  this  case  |Pn+i|n+i|  behaves  like  aqir 22^aq t2 ,  i.e  it  is  increasing  for  large  t.  The  function  has 
a  pole  at  t  =  0  and  for  t  >  0  it  is  positive.  So  that  it  has  a  local  minimum  for  t  >  0.  The 
question  ism  does  it  have  a  local  minimum  at  t  >  1.  This  largely  depends  on  a. 
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An  example  is  shown  in  Figure  16.  In  this  example  q  =  1,  R  = 


Figure  16:  Determinant  of  State  Covariance  as  a  function  of  scheduling  time  when  probability 
of  detection  is  improved  by  scheduling. 


6.3.3  Size  and  Shape 

•  Scheduling  improves  range  measurement: 

(30) 


(31) 

Further  calculations  show,  that  in  this  situation  |Pn+1|ri+1|  is  decreasing,  has  no  local 
minima  in  the  domain,  and  behaves  like  c\C2/t  for  large  t.  This  means  that  scheduling 
time  always  results  in  information  gain.  Further  calculations  show,  that  in  this  situation 
|-Pn_l_i|n_l-i|  is  decreasing,  has  no  local  minima  in  the  domain,  and  behaves  like  cic2/t  for 
large  t.  This  means  that  scheduling  time  always  results  in  information  gain. 

Numerical  examples  are  represented  in  Figure  17. 

6.3.4  Orientation  of  the  Measurement  Noise  Covariance  Changes  with  Scheduling 
Time 

Kalman  covariance  update  can  be  re-written  as 

Pn+l\n+l  =  (Pn+l\n  +  R  *)  ^  (32) 

Suppose  it  is  possible  to  construct  a  waveform  so  that  R  in  the  form  R  =  ARqAt  ,  where  Ro  is 
fixed  and  diagonal  with  the  elements  on  the  diagonal  arranged  is  in  descending  order,  and  A  is 


R  = 


C2&2  +  Cl/f,  bC2 


bC2 


C2 


Scheduling  improves  Doppler  measurement: 

( (c262)/t  +  ci  ( bc2)/t 

V  {bc2)/t  c2/t 
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Scheduling  Time 

Figure  17:  Determinant  of  State  Covariance  as  a  function  of  scheduling  time  when  only  range 
or  doppler  measurement  is  improved.  In  this  example  q  =  1,  b  =  1,  ci  =  1,  C2  =  1 

a  rotation  matrix  with  known  rotation  angle.  Fielder  [20]  tells  us  that  [Pn_)_1|n+1|  is  maximized 
iff  Pn+i\n  =  AEAT ,  where  E  is  a  diagonal  matrix  with  element  on  the  diagonal  (eigenvalues) 
arranged  in  the  ascending  order.  In  2D  this  makes  covariance  ellipses  perpendicular  to  each 
other.  Lets  assume  that  it  takes  time  T  >  1  to  generate  R  =  R*  perfectly  perpendicular  to 
Pn+i\n-  The  question  is:  if  we  use  R*  at  times  t  <T,  do  we  improve  the  accuracy  of  tracking  by 
taking  a  measurement  at  this  time.  Lets  first  look  at  an  example.  Lets  say  that  at  time  T  =  10 
the  ellipse  of  predicted  covariance  matrix  Pn+i\n  main  axis  was  0  degrees  with  x  -axis  (range), 
then  the  waveform  is  constructed  so  that  its  covariance  matrix  is 

R  =  A(4>)  ^ A(0)_1,  n  <  r2,  (33) 

then  |Pn+1|n+1|  is  the  smallest  possible  with  respect  to  4>.  Is  it  also  the  smallest  with  respect 
to  f?  Numerical  simulations  show,  that  the  answer  to  the  question  is  negative.  This  happens 
because  the  gain  caused  by  orientation  optimization  is  negligible  compared  to  the  loss  caused 
by  the  growth  of  Pn+\\n  with  time.  An  example  is  shown  in  Figure  18.  This  example  suggest 
that  measuring  at  time  t  =  1  is  best,  even  though  the  measurement  covariance  matrix  is  “under¬ 
rotated”  . 

6.4  Conclusion 

We  have  investigated  several  situations  to  answer  the  question:  how  the  time  spent  on  scheduling 
the  radar  resources  (to  achieve  better  measurement)  affects  the  accuracy  of  the  target  state. 
We  postulated  that  the  measurement  quality  improves  by  spending  longer  time  at  scheduling. 
The  longer  time  between  measurements,  on  the  other  hand,  causes  larger  spreading  in  target 
uncertainty.  We  have  computed  the  quantitative  values  of  the  parameter,  for  which  scheduling 
improves  target  estimation  and  gains  more  information. 

We  also  have  described  situations  in  which  no  improvement  is  possible  and,  in  fact,  in  these 
situations  scheduling  will  result  in  loss  of  accuracy  and  consequently  information. 
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3.5 


Figure  18:  Determinant  of  State  Covariance  as  a  function  of  scheduling  time  when  the  orientation 
of  the  measurement  covariance  is  improved.  In  this  case  scheduling  will  cause  the  loss  of  accuracy 

7  More  Simulation  results 

Recently,  we  have  further  investigated  two  problems,  i.e.,  waveform  scheduling  and  Doppler 
resilient  Golay-type  pulse  trains  and  these  works  have  produce  close  to  submission  papers.  In 
addition,  we  also  done  some  experimental  simulations  in  line  with  the  waveform  scheduling. 

7.1  Waveform  Scheduling 

Adaptive  scheduling  of  radar  and  sonar  sensor  resources  has  been  shown  by  many  studies  to 
produce  significant,  and  some  times  enormous,  gains  in  performance  [21,  22].  Generally,  the 
best  performance  is  achieved  by  scheduling  over  an  extended  horizon  since  this  allows  the  future 
consequences  of  current  actions  to  be  accounted  for.  However,  finite  horizon  scheduling  is 
computationally  expensive  and  thus  can  rarely  be  used  in  practice.  In  this  work  we  have,  in 
collaboration  with  researchers  at  Colorado  State  University,  specifically,  Drs  Ali  Pezeshki  and 
Edwin  Chong  and  their  student  Zhenliang  Zhang,  sought  ways  of  guaranteeing  the  performance 
of  computationally  efficient,  but  sub-optimal,  greedy  scheduling  schemes.  In  particular,  we 
have  made  significant  advances  in  the  area  of  finite  horizon  scheduling  with  string  submodular 
objective  functions.  This  will  enable  performance  bounds  on  longer  term  scheduling  of  waveforms 
for  radar,  and,  in  particular,  hopefully  show  that  greedy  methods  for  scheduling  waveforms  can 
provide  close  to  optimal  performance. 

Previously  we  showed  that  for  objective  functions  which  are  string  submodular  a  greedy 
approach  to  choosing  a  sequence  of  actions  achieves  at  least  a  (1  —  e_1  ^approximation  compared 
to  the  optimal  approach.  More  recently,  we  have  shown  that  these  bounds  can  be  improved  by 
considering  further  properties  of  the  objective  function.  In  particular,  if  the  objective  function 
has  total  backward  curvature  a  than  the  performance  bound  for  greedy  scheduling  becomes 
(1  —  e~a)/a  of  the  optimal  strategy  while  for  total  forward  curvature  e  the  bound  is  (1  —  e).  These 
results,  which  will  appear  in  a  forthcoming  paper,  more  tightly  bound  the  loss  in  performance 
for  using  computationally  efficient  greedy  scheduling  in  place  of  optimal  scheduling  [23]. 

We  are  currently  investigating  the  application  of  these  general  results  to  particular  problems 
in  radar  and  sonar  tracking. 
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7.2  Doppler  Resilient  Golay-type  Pulse  Trains 

Suppression  of  sidelobes  in  the  ambiguity  function  is  essential  for  detection  of  closely-separated 
objects  in  radar  and  sonar  processing.  It  is  well-known  that  complementary  pulse  trains  are 
useful  for  sidelobe  suppression.  Previously,  we  have  shown  that,  to  achieve  this  suppression  at 
non-zero  Dopplers,  the  sequences  must  be  sent  in  a  particular  order.  With  Dr  Ali  Pezeshki  and 
his  student  Wenbing  Dang,  we  have  recently  extended  these  results  by  shaping  the  pulse  trains 
at  both  the  transmit  and  receive  ends.  Using  a  receive  pulse  train  which  is  not  matched  to  the 
transmit  pulse  train  allows  an  extra  degree-of-freedom  in  shaping  the  ambiguity  function,  albeit 
with  some  loss  in  SNR. 

The  benefits  of  independently  scheduling  transmit  and  receive  pulse  trains  are  shown  in 
Figure  19.  The  scene  has  three  strong  reflectors  and  two  weak  reflectors.  The  conventional  ap¬ 
proach,  in  which  alternating  Golay  complementary  waveforms  are  used  for  transmit  and  receive, 
performs  poorly  with  the  weak  reflectors  suppressed  by  the  strong  reflectors.  Performance  is 
significantly  improved  using  the  proposed  binomial  transmit-receive  pulse  sequence  design.  It 
can  be  seen  that  the  binomial  design  provides  sidelobe  suppression  for  Dopplers  in  the  range 
[—1,1]  rad.  Compared  to  existing  techniques  [24],  this  increases  by  a  factor  of  10  the  range 
of  sidelobe-suppressed  Dopplers.  As  a  result,  weak  reflectors  are  not  masked  by  the  strong 
reflectors  in  this  example.  These  results  are  contained  in  the  forthcoming  paper  [25]. 


<  10~5  N=16,  Conventional  design,  output  in  dB 


N=16,  Binomial  design,  output  in  dB 
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Figure  19:  Delay-Doppler  gain  using  (a)  conventional  and  (b)  binomial  transmit-receive  design. 


7.3  Complementary  Pair  Scheduling 

Complementary  waveforms  achieve  lowest  sidelobes  at  zero  Doppler.  The  idea  is  to  sched¬ 
ule  transmission  order  of  waveforms  to  achieve  optimal  off-zero  Doppler  performance.  This  is 
achieved  by  using  a  Reed-Miiller  sequence  tuned  to  Doppler  region  to  provide  transmission  or¬ 
der.  We  have  carried  out  simulations  based  on  our  recent  work  in  [17]  and  results  are  shown 
in  figures  20  and  21,  which  demonstrated  that  using  a  scheduled  Reed-Miiller  sequence,  20dB 
improvement  in  SNR  and  30%  improvement  in  probability  of  detection  are  gained  over  the 
standard  LFMs  and  standard  complementary  sequences. 

7.4  Beam  scheduling 

We  also  simulated  the  situation  of  a  steerable  phased  array  radar.  The  aim  is  to  selected  the 
pointing  direction  so  as  to  combine  tracking  and  search.  We  use  the  cardinalised  probability 
hypothesis  density  (CPHD)  for  filtering  and  attempt  to  maximise  the  probability  of  target 
detection.  The  CPHD  maintains  tracks  as  well  as  potential  targets  in  unobserved  areas.  The 
detection  probability  is  an  inexpensive  proxy  for  information  gain.  As  shown  in  Fig.  22,  the 
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(a) 


(b) 


Figure  20:  (a)  Sidelobes  for  standard  LFMs,  standard  complementary,  and  Reed-Muller  se¬ 
quences.  (b)  SNR  for  standard  LFMs,  standard  complementary,  and  Reed-Muller  sequences. 


Time  (sec) 

Figure  21:  Detection  probability  for  standard  LFMs,  standard  complementary,  and  Reed-Muller 
sequences. 


time  required  for  measurement  collection  can  be  reduced  by  40%  by  using  this  beam  scheduling 
compared  to  regular  surveillance. 

7.5  LFM  waveform  scheduling 

Another  piece  of  work  we  performed  is  the  simulation  analysis  of  the  LFM  waveform  scheduling. 
It  is  understand  that  scheduling  over  long  horizon  is  computationally  expensive.  On  the  other 
hand,  if  tracking  is  performed  with  a  filter  of  Kalman  type,  the  long  horizon  scheduling  is 
performed  off-line.  The  objective  function  to  be  designed  is  such  that  selects  waveforms  from 
a  library  to  maximize  the  information  gain.  Clearly,  the  optimisation  policy  depends  on  the 
dynamics  of  the  underlying  target.  In  the  simulation,  two  examples  of  dynamics  were  considered: 
1  -  random  walk  with  added  range  uncertainty  (i.e  sea  targets),  2  —  random  walk.  The 

simulation  result,  shown  in  Fig.  23,  pre-computing  the  optimal  sequence  of  LFM  waveforms  has 
1/m^th  the  expense  of  scheduling  over  a  horizon  h  with  m  candidate  waveforms. 
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Figure  22:  (a)  Revisit  counts  for  for  regular  and  schedued  surveillance,  (b)  RMSE  for  regular 
and  schedued  surveillance. 


Figure  23:  Optimal  policies  for  two  different  process  noise  covariance  matrices. 
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